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[Article  by  0.  K.  Voronkov  and  V.  I.  Marov] 

1.  The  Task  of  the  Investigations 


[Text]  To  calculate  the  dynamic  elastic  modulus  E , Poisson's  ratio  /u  , the 
shear  modulus  G,  and  ocher  characteristics  of  the  ^elasticity  of  rocks'^it  is 
necessary,  besiues  the  longitudinal  wave  velocity  v , to  know  the  transverse 
wave  velocity  vs  or  the  Rayleigh  wave  velocity  v...  p Whereas  the  determin- 
ation of  vg  and  v^  during  seismoacoustic  observations  in  mountain  workings 
(shafts  and  wells)  causes  no  difficulties  in  principle,  during  engineering 
seismic  surveying  on  the  day  surface  the  determination  of  vg  and  vR  involves 
a number  of  difficulties.  In  particular,  in  a permafrost  region  the  pre- 
sence of  a high-velocity  l^yer  of  frozen  unconsolidated  deposits  in  situ  on 
bedrock  excludes  the  possibility  of  obtaining  P^^P^  and  waves  (ex- 

change on  the  boundary  of  unconsolidated  and  solid  rocks)  — intensive  oscil- 
lations, registered  with  certainty  in  a number  of  regions  where  there  is  an 
absence  of  frost  [ 7 a 3] - It  is  impossible  to  form  those  types  of  waves  at 


The  nature  of  waves  registered  in  each  new  region  of  investigations  is  deter- 
mined by  a large  group  of  kinematic  and  acoustic  characteristics,  namely: 
the  absolute  wave  velocities,  the  ratios  of  the  velocities  of  different  waves, 
the  relative  intensity,  the  ratios  of  the  oscillation  periods,  the  absorption 
coefficients  and  decrements,  etc.  By  constructing  variational  curves  or  histo- 
grams for  a statistical  series  of  each  of  those  signs  one  can  determine  their 
average  values,  and  also  the  basic  range  of  variation  of  the  corresponding 
characteristics.  In  addition,  other  characteristics  of  the  recording  of  waves 
are  used:  the  traceability  in  the  region  of  first  or  subsequent  arrivals,  the 

characteristic  shape  of  the  travel-time  curve,  etc.  Comparisons  with  theo- 
retical seismograms,  etc,  are  also  used.  The  nature  of  the  waves  is  usually 
determined  for  the  region  of  the  work  as  a whole.  As  a rule,  the  nature  of 
the  waves  cannot  be  established  for  local  parts  of  that  region.  Therefore 
it  is  very  important  to  know  how  given  geological  factors  which  vary  within 
a given  region  (the  thickness  of  deposits,  their  thawed  or  frozen  state,  etc) 
influence  the  properties  of  the  wave  field. 
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Seismic  prospecting  investigations  (engineering  and  ore  seismic  prospecting) 
in  a permafrost  region  have  established  that  the  most  intense  wave  on  the 
seismograms  is  the  Rayleigh  wave  registered  in  "second  arrivals"  — after 
longitude  refracted  waves.  With  respect  to  an  overwhelming  majority  of  the 
cases  of  engineering  and  ore  seismic  prospecting  it  can  be  regarded  as  a 
Rayleigh  wave  in  a layer  (unconsolidated  deposits)  on  a half-space  (bedrock 
is  solid  rock).  The  phase  velocity  of  such  a wave  depends  on  a number  of 
factors  (the  elasticity  and  density  of  the  layer  and  half-space,  the  layer 
thickness  and  the  wave  length.  Therefore,  before  using  the  experimental 
values  of  v^  obtained  in  field  seismic  prospecting  investigations,  for  the 
calculation  of  E^,  ^ , etc,  of  unconsolidated  and  solid  rocks  it  is  necessary 
to  take  into  consideration  quantitatively  the  influence  of  the  enumerated 
factors  on  the  phase  velocity  v^. 

In  order  to  study  the  wave  field  registered  during  the  seismic  prospecting 
of  small  depths  in  regions  of  the  Extreme  North  and  permafrost  an  ultrasonic 
modeling  was  carried  out  on  a three-dimensional  model.  Besides  that,  to 
study  the  dispersion  of  the  velocity  of  Rayleigh  waves  in  a layer  on  a half- 
space the  theoretical  dispersion  curves  were  calculated  and  an  analysis  was 
made  of  the  field  seismic  prospecting  material  obtained  in  permafrost  r ;ions 
of  the  Northeast  of  the  USSR  (the  basins  of  the  Indigirka,  Vilyuy,  Kolyma 
and  Tatta  rivers). 

2.  Characteristics  of  the  Wave  Picture  According  to  Data  of  Three-Dimensional 
Ultrasonic  Modeling  (Layer  on  a Half-space) 

To  study  the  wave  picture  with  reference  to  the  conditions  of  permafrost  re- 
gions the  authors  conducted  an  ultrasonic  modeling  of  the  wave  picture  by  the 
method  of  longitudinal  profiling  zz.  The  profiling  was  carried  out  both  from 
the  surface  of  the  layer  and  from  the  surface  of  the  half- space  (Figure  l). 

In  the  general  case  the  three-dimensional  model  represented  a layer  (or  two 
layers)  of  different  thickness,  composition  and  state  on  the  half-space.  The 
materials  for  the  making  of  models  were  natural,  particularly  during  the  cre- 
ation of  layers  use  of  made  of  water  (v^,  = 1450  m/sec,  density  o = 1 g/cc), 
ice  (vp  = 3650  m/sec,  <5  = 0.91  g/cc),  frozen  strongly  icy  loam  (v?  = 3700 
m/sec,  = 2000  m/sec,  cf  = 1.9  g/cc,  and  frozen  ice-saturated  ssmd  (v  = 

= 4000  m/sec,  v^  = 2400  m/sec,  6 = I.85  g/cc).  P 

The  half-space  is  a granite  block  set  on  porolon  and  placed  in  a hydrolyzed 
vessel.  The  block  dimensions  were  430  x 190  x 80  mm;  the  wave  velocities  along 
the  face  length  were:  Vp  v^ 

- in  the  dry  state  4210  m/ sec  2650  m/ sec 

in  the  water-saturated  state  5350  " 2820  " 

- in  the  ice-saturated  state  6000  " 3200  " 

Let  us  examine  the  questions  of  similarity  of  the  model  (m)  and  nature  (n). 

VP(H>  VS  * M * 

As  is  known  [5],  the  velocity  criterion  of  similarity  is  C c 

which  means  equality  of  the  Poisson's  factors  = )*(my  density  cri- 

f.\l, 

terion  • — : . For  the  models  in  our  experiments  Cv  = 1 and  C<f  = 1 and, 
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Figure  1.  Principal  diagrams  of  three-dimensional  ultrasonic 
modeling  of  elastic  waves  with  a layer  (layers)  on  a half- 
space: 1 - ice;  2 - granite;  3 - loam;  4 - sand;  5 - water; 

I — radiator;  R — receiver  of  ultrasonics  a - Frost;  b - 
Thawed  water  on  frost. 


according  to  B.  N.  Ivakin,  such  models  belong  in  the  first  group  and  for  them 
the  constants  of  similarity  of  the  frequency  time  depends  only  on  the  constant 
of  geometric  similarity: 


We  are  studying  a layer  on  a half-space . A body,  the  linear  dimensions  of 
which  are  greater  than  (2-3)  A,  where  A is  the  wave  length,  can  be  considered 
a half-space  for  the  study  of  seismic  processes.  For  purposes  of  engineering 
seismic  prospecting  in  nature,  of  interest  is  the  case  of  variation  of  the 
layer  thickness  K = 0 - 100  m (h/A  = 0-3).  If  in  the  model  H = 0 - 10  cm 
(H/A  = 0-3),  then  C = 1000;  Ca  = 10-3.  Tf  in  nature  ^ ^ 100  Hz,  then  in 
work  on  the  models  tne  wave  frequencies  must  be  of  the  order  of  100  kHz. 

During  ultrasonic  modeling  an  IPA-59  seismoscope  was  used,  and  the  radiators 
and  receivers  (hydraulically  and  thermally  insulated)  had  a natural  frequency 
of  150-200  kHz.  The  visible  frequencies  of  the  waves  on  the  recording  varied: 
100-130  kHz  (Rayleigh  and  transverse  waves  in  granite  and  frozen  layers)  and 
170-180  kHz  (a  forward  wave  in  water).  Experiments  on  frozen  models  were  con- 
ducted in  winter  in  the  open  air  at  from  -10  to  -15°C.  For  ultrasonic  profil- 
ing the  radiator  was  set  at  the  zero  point  and  the  receiver  was  moved  with  a 
step  of  2 cm.  The  radiator  and  receiver  were  moistened  with  water  and  frozen 
to  the  model,  and  that  assured  good  acoustic  contact.  Two  recordings  were 
obtained  at  each  point:  with  a large  amplification  of  the  registering  appar- 

atus (distinct  first  arrivals  of  waves)  and  lesser  amplifications  (resolved 
recording  of  subsequent  waves).  The  obtaining  of  seismograms  with  "strong" 
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Figure  2.  a - Oscillogram  of  ultrasonic  profiling  on  the 
surface  of  frozen  sand,  H =4  cm;  b - Oscillogram  of  ultra- 
sonic profiling  on  the  surface  of  ice,  Hice  = 5 cm. 
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Figure  3.  Oscillogram 
of  ultrasonic  profiling 
on  the  surface  of  water, 
\ = 2 ci. 


and  "weak"  recordings  is  characteristic  of  the 
method  of  engineering  seismic  prospecting.  The 
wave  picture  was  photographed  from  the  screen  of 
the  seismoscope.  Examples  of  recordings  are 
presented  on  Figures  2 and  3.  Presented  below  is 
a description  of  the  main  characteristics  of  the 
wave  picture  according  to  the  data  of  modeling  for 
conditions  characteristic  of  seismic  prospecting 
in  regions  of  the  Extreme  North  and  permafrost 
(Figure  l).  Examples  of  travel  time  curves  of  the 
principal  types  of  waves  are  shown  on  Figures  4- 5. 


Figure  4.  Travel  time  curves  of  the  principal 
types  of  waves  in  profiling  zz  on  frost.  1-3*. 
profiling  on  a half- space  loaded  ’with  a layer 
(diagram  b.  Figure  l);  4 - profiling  on  a half- 
space of  loam  (diagram  c.  Figure  l);  5 - profil- 
ing on  the  surface  of  sand  (diagram  d,  Figure  l). 
a - microseconds 


A layer  of  ice  on  a frozen  half- space  (profiling 
on  the  surface  of  ice,  Figure  la).  At  a layer 
thickness  H = 0-1  cm  the  recordings  are  strong 
and  not  complicated  by  interference : 2 or  3 period 
(2-3 )T  oscillation  of  the  longitudinal  wave 
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in  granite  (v  = 5900-6100  m/sec,  ^ = 125-150  kHz)  and  (2-3)  T of  oscil- 
P2 

lation  of  the  Rayleigh  R wave  (v  = 3200  m/sec  at  H = 0 cm,  v„  = 3000  m/sec 

at  H = 0.5  cm,  v^  = 2950  m/sec  at  H = 1 cm;  £ = 100-120  kHz).  The  ratios 
of  the  maximal  amplitudes  of  oscillations  are’: 

a)  /A  = 10-20  (at  H = 0 cm)  in  both  the  air-dry  and  the  ice-saturated 
^ P2 

blocks: 

b ) A_/A  = 5 '8  ( at  H = 0.5  cm) ; 

R ?2 

c)  A/A  = 2-3  (at  H = 1 cm). 

R Po 


Thus  the  presence  of  a layer  of  ice  even  of  relatively  small  thickness 
H/A^ir  l/2  greatly  reduces  the  amplitude  of  Rayleigh  waves  and  to  a lesser 
degree  vR  (by  6-8%). 


Figure  5.  Travel  time  curves  of  the  principal  types  of  waves 
during  profiling  zz  on  frost,  on  the  surface  of  a layer  (dia- 
gram of  Figure  la).  a - microseconds 
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At  2 era  on  recordings  of  profiling  are  visible  (l-2)  T oscillations  of 
a forward  longitudinal  wave  in  ice  (v  = 3500-3700  ra/sec,  ■/  = 170-190  kHz), 
a longitudinal  refracted  wave  in  P1  granite  (v  = 6COO-6lOO  ra/sec,  = 

= 125-160  kHz)  and  a Rayleigh  wave  in  a layer  of  '2  ice  (vR  = 1700-1800 
ra/sec,  f = 100-125  kHz),  the  amplitude  of  which  is  10-15  "1  times  as 

large  as  A^  . An  interesting  feature  of  the  oscillograms  is  the  presence  on 
many  of  ‘2  them  of  recordings  of  multiple  reflected-refracted  waves  (1-2)  T 
of  oscillations  with  v = 5900-6100  ra/sec,  /■  = 125-140  kHz,  the  amplitude  of 
which  is  2-4  times  greater  than  A . The  wave  was  registered  in  different 
arrangements  of  observations  p2  (Figure  l)  at  distances  of  more  than  6-10 
cm  from  the  radiator,  with  receivers  satisfying  single  or  double  reflection 
from  a free  surface. 


An  ice  layer  on  a frozen  half-space  (profiling  on  the  surface  of  the  half- 
space). Three  principal  waves  were  found  on  the  recordings:  longitudinal  in 
granite  --  (l-2)  T oscillations  with  v = 5900-6100  ra/sec,  -j~  - 125-140  kHz; 
a longitudinal  wave  in  ice  — 2T  p2  oscillations  with  Vp  = 3700  ra/ sec , 

the  amplitude  of  which  is  1.5-2  times  larger  than  A ; a 1 Rayleigh 
wave  with  ■f'  = 110-125  kHz.  Ap/A^  = 5-10  at  p2  distances  near  the 

radiator  and  A /A^  = 2-5  at  *2  'distances  remote  fr ora  the  radiator.  The 

velocity  vR  p rt2  = 3000  ra/sec  (Hal  cm)  and  on  travel  time  curves  of 
different  ‘phases  of  the  R wave  vR  = 2100-2U00  ra/sec  (H  = 10  cm). 


Thus  in  the  presence  of  a layer  on  a half- space  a Rayleigh  wave  registered  on 
the  surface  of  the  half-space  is  damped  with  distance  considerably  more  than 
on  a free  half-space,  and  vR  decreases  with  increase  of  H. 

A layer  of  frozen  sand  and  loam  on  a frozen  half-space  (profiling  on  the  sur- 
face of  the  layer).  The  wave  picture  with  respect  to  frequencies  and  ampli- 
tudes is  similar  to  the  above-considered  case  of  an  "ice  layer  on  a half -space. 
However,  during  profiling  on  sand  and  loam  stronger  damping  of  the  field  with 
distance  is  observed  and  a recording  of  oscillations  less  stable  in  form  in 
comparison  with  the  profiling  on  ice.  Visible  on  the  recordings  are  (2-3)  T 
oscillations  of  a forward  wave  in  the  layer  (sand  — vp  = 3900-4050  ra/sec, 
loam  — Vp  = 3700  m/sec,  -f  = 100-140  kHz),  of  a refracted  wave  in  granite 

with  vp2  = 6OOO-6IOC  m/sec,  -f  = 125-140  kHz.  A /A  2-3;  of  a Rayleigh 
wave  in  a bed  (sand  — vR  = 2400  ra/sec,  loam  P1  p2  vR  = 1950-2000  m/sec, 

A_  /A  = 10-20)  and,  starting  with  24-26  cm  from  the  radiator,  a Rayleigh 
*1  p2 

wave  with  vR  = 2950  m/sec  and  an  amplitude  2-5  times  greater  than  A^  . At 
greater  distances  from  the  radiator  there  is  no  deterioration  cf  '2  the 
Rayleigh  wave  packet  (the  phase  velocity  of  different  extremums  corresponds  to 
vR  of  sand  or  loam). 


A water  layer  on  frost  (Figure  1,  e and  f).  A characteristic  feature  of  the 
wave  picture  is  a predominance  of  longitudinal  waves:  in  ice  with  v = 3600- 
3700  ra/sec  (at  negative  temperatures)  and  v = 3200-3300  m/sec  ^(melting 
ice),  j"  = 140-200  kHz;  in  granite,  refracted  and  reflected-refracted  with 
v = 5800-6100  m/sec,  f = 125-140  kHz  and  in  water  with  vp  = 1450-1500  m/sec, 
p2  p 
f = 100-110  kHz. 
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Whereas  the  first  two  waves  are  of  approximately  identical  intensity,  the 
amplitude  of  that  in  water  is  approximately  8-10  times  greater.  In  a layer 
of  ice  of  5 cm  and  during  profiling  on  the  surface  of  the  water  an  exchange 
PSP  wave  also  is  positively  registered,  propagating  over  the  ice  with  a phase 
velocity  of  1900  m/sec,  j'  = 110-125  kHz. 

A thawed  layer  on  a thawed  half-space  (Figure  1,  g,h,i).  In  water- saturated 
granite  Vp  = 5300  m/sec  and  v = 2810  m/sec.  During  profiling  on  the  surface 
of  water  Figure  lg)  (2-3)T  oscillations  were  registered;,  a very  intense  for- 
ward longitudinal  wave  in  water  (v  = 1^50-l^CO  m/sec,  ■f'=  100  kHz,  refracted 
and  reflected-refracted  waves  corresponding  to  the  surface  of  granite  (v  = 

= 5350-5^00  m/sec,  j~  = 120-lhO  kHz)  and  an  exchange  wave  PSP  ~2 

v = 2900  m/sec,  f = 100-120  kHz.  The  amplitude  of  the  latter  is 

6-12  times  A , The  amplitude  of  the  reflected-refracted  wave  exceeded  A_ 

P2  P ?2 

t>y  2-3  times.  During  profiling  on  the  surface  of  granite  (Figure  Ih)  (2-3)  T 
oscillations  were  registered:  a longitudinal  (v^  = 5200-5^-00  m/s ec,  -j-  = 125- 

150  kHz),  a reflected-refracted  with  a greater  “2  amplitude  than  Ap  , and  also 
a Rayleigh  wave.  The  phase  velocity  of  the  R wave  diminished  ' 2 slightly 

with  increase  of  H:  from  2810  m/sec  (H  = 0 cm)  to  2690  m/sec  (H  = 10  cm). 

The  loading  of  the  half-space  with  a layer  of  water  has  a great  influence  on 
the  change  of  A^:  whereas  at  H = 0 AR/Ap^  = 8-10,  at  H = 5-10  A^/A^  = 2-k, 

and  in  the  presence  of  the  layer  the  damping  of  the  Rayleigh  wave  increases 
with  distance. 


During  profiling  on  the  surface  of  granite  in  the  presence  of  a layer  of 
water-saturated  sand  on  it  (Figure  li),  as  in  the  case  considered  above,  a 
longitudinal  and  a Rayleigh  wave  were  registered  in  granite.  With  increase 
of  the  layer  thickness  V_  decreased:  2810  m/sec  (H  = 0 cm),  2700  m/sec  (H  = 

= 1 cm),  2330  m/sec  (H  = 5 cm)  and  2100  m/sec  (H  = 10  cm).  The  ratio  A_/A 
during  the  loading  of  the  half-space  with  a sand  layer  decreased  to  2-5^  ^ 

(at  H * 5-10  cm)  and  the  R-wave  damping  increased  with  distance. 


3.  Study  of  the  Dispersion  of  the  Velocity  of  Rayleigh  Waves  in  a Layer 

on  a Half-space  With  Reference  to  Problems  and  Conditions  of  engineering 
Seismic  Prospecting 


Calculating  Dispersion  Curves 

We  will  examine  a layer  (medium  l)  on  a half-space  (medium  2).  It  is  known 
from  seismology  [h,  etc!  that  the  phase  velocities  of  short  Rayleigh  waves 
approach  v in  a layer  (since  short  waves  do  not  penetrate  deeply),  and  the 
phase  velocities  of  very  long  waves  approach  v„  of  the  half-space.  Let  us 
turn  to  a solution,  known  from  the  theory  of  elasticity  [6],  of  the  problem 
of  propagation  of  transverse  waves  in  a layer  with  the  thickness  H on  a half- 
space: Q 

7(V/)  = n7, 
ivi 


where 


*\  S1 
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are  the  shear 


uj  = P.  njl  is  the  circular  frequency;  T is  the  period;  G.  and 

moduli,  v and  v are  the  transverse  wave  velocities  in  the  layer  and  the 
S1  S2 

half-space;  v is  the  transverse  wave  velocity  in  the  layer  on  the  half-space 
(along  the  surface  of  the  layer). 


Equation  (l)  determines  in  implicit  form  the  dependence  of  v on  the  frequency; 
it  has  solutions  only  at  real  values  of  and  mo,  so  that  “it  always  happens 

that  v >v  >v  . Therefore  the  propagation  of  ohe  type  of  wave  under  consii- 
s2  s 

eration  is  oossible  only  at  v > v . We  will  write  (l)  in  Jne  expanded  form: 

s2  S1 


tV* 


where 


s 


i 


are  the  densities  of  the  layer 


and  the  half-space. 


All  values 


, where  n = 1,2,  ...,  are  discarded. 


For  the  conditions  of  seismic  prospecting  of  small  depths  the  ratio  v^/v„ 
varies  in  a relatively  narrow  range  (0.91-0.9*0,  and  under  permafrost' 
conditions  (a  frozen  icy  layer  on  the  same  half-space)  in  still  narrower  limits 
(0.925-0.935).  Therefore  we  will  write  (la)  in  a form  convenient  for  further 
analysis  and  comparison  with  experiment: 


" *>" 
VK, 


Equation  (2),  determination  of  v , was  solved  on  an  M-20  computer  with  the 
following  values  of  the  parameters  included  in  it: 


- 1 .5;  v^  = 1200  m/sec  - 4000  m/sec  (steps  of  400  m/sec); 

v = 100  - 2000  m/sec  (steps  of  200  m/sec)  for  the  values  vn  >v_; 
R & 

H/A  = 0 - 10  (at  H/A  = 0 - 0.5  a step  of  0.05,  at  H/\  = 0.5  - 1 a 

step  of  0.1  and  at  H/A  = 1 - 10  a step  of  l). 


Examples  of  dispersion  curves  are  shown  on  Figures  6-8.  At  the  relatively  low 
values  = 0.1  - 0.5  km/sec  (which  is  characteristic  of  thawed  unconsolidated 
deposits1^-)  the  values  of  v decrease  sharply  from  vn  to  vD  during  change  of 


H/A  from  0 to  0.2-0. 3.  Characteristic  of  a frozen  layer  (v  = 1. 5-2.0  km/sec) 
on  a half -space  is  a smoother  change  of  v from  vn  to  vw  K1  at  H/A  =0-1. 
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Figure  6.  Dispersion  of  Rayleigh  wave 
velocities  v^  in  the  layer  on  a half- 
space.  vR  — the  velocity  of  Rayleigh 
waves  in  Ithe  layer,  km/sec;  v — 
the  velocity  of  Rayleigh  waves  K2 
in  the  half-space,  km/sec. 


Figure  7.  Theoretical  dispersion 
curves  of  the  Rayleigh  wave  velocity 

v in  the  layer  on  a half-soace. 

R 

v and  v as  on  Figure  6. 

R1  R2 


Experimental  Data  on  the  Dispersion  of  the  Velocity  of  Rayleigh  Waves 
in  Permafrost  Regions 

Rayleigh  waves  registered  during  shallow  seismic  prospecting  in  a permafrost 
region  have  been  described  in  a number  of  works  [l,  2].  Characteristic  of 
each  region  is  a relatively  narrow  range  of  their  frequencies  and  consider- 
able changes  of  frequency  of  R-waves  are  observed  only  in  the  transition  of 
a wave  from  a region  of  frosted  rocks  into  a region  of  thawed  rocks  [2],  which 
is  a reliable  sign  for  the  discovery  of  taliks.  The  amplitude  of  R-waves,  as 
a rule,  are  5-15  times  as  large  as  the  amplitudes  of  forward  and  refracted 
longitudinal  waves.  Presented  below  are  seismic  prospecting  data  on  the 
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Figure  8.  Dispersion  of  the  velocity  of 
Rayleigh  waves  v„  in  the  presence  of  a 
layer  on  a half-space:  a - calculating 
dispersion  curves  at  different  v in 
the  hall-space  and  at  a v_,  of  A2 
the  layer  of  1.?  km/sec;  1 b - disper- 
sion of  v^  in  a frozen  layer  cn  a frozen 
half-space:  1 - theoretical  curve  (ice 
on  granite);  2 - experimental  curve  (ice 
on  granite)  according  to  ultrasonic  model- 
ing; 3,  4 & 5 - data  of  seismic  prospect- 
ing in  permafrost  regions:  basins  of  the 
rivers  Vilyuy,  Indigirka  and  Kolyma  re- 
spectively; 6 - dispersion  of  vR  in  the 
semi-space  (ultrasonics  on  the  surface  of 
granite)  in  the  presence  of  layers  of 
water,  wet  sand  and  ice  on  it. 


dispersion  of  v^  in  permafrost  regions  (Table  l).  The  values  of  v^  were  cal- 
culated by  the  method  of  the  difference  travel  time  curve,  which  made  it  pos- 
sible to  exclude  the  influence  on  v^  of  difference  in  the  installation  of 
seismic  receivers  and  local  heterogeneities  in  the  upper  part  of  the  profile. 
To  exclude  the  influence  on  v^  of  heterogeneities  of  a different  type  in  the 
layer  and  half-space  (variations  of  granulometeric  composition,  iciness,  tem- 
perature in  the  layer,  and  differences  in  the  degree  of  fracturing  of  the 
half-space  — the  bedrocks)  within  the  region  of  investigations  in  the  given 
range  of.  value  of  H we  calculated  the  average  weighted  value  of  v^  along  the 
survey  lines.  Due  to  difference  in  the  length  of  the  survey  lines  for  each 
range  of  H values  the  calculated  values  of  vR,  generally  speaking,  are  not  of 
equal  accuracy.  The  values  of  v were  determined  on  sections  with  minimal 
values  of  H/^.  As  a rule  v^  is^5-10%  lower  than  v^  in  preserved  ice-satur- 
ated samples  of  corresponding  rock  at  the  given  temperature. 


The  Vilyuyskaya  GES.  The  base  was  composed  of  fine-grained  and  medium-grained 
dolerites  of  the  Lower  Mesozoic  (T^).  The  thickness  of  the  belt  of  weathering 
of  the  recks  is  from  0 to  30  meters.  Unconsolidated  thick  deposits,  as  a rule, 
consist  in  the  first  few  meters  of  loams,  more  rarely  of  sand  loams,  clays  and 
sands.  The  thickness  of  the  permafrozen  rocks  in  the  banks  of  the  Vilyuy  River 
valley  is  150-200  meters.  The  temperature  t of  the  permafrozen  rocks  on  the 
investigated  section  at  depths  of  15-20  meters  is  about  -2  to  -3°C.  The  aver- 
age characteristics  of  the  half-space  (the  rocks  below  the  belt  of  weathering) 
are : v = 5000  m/sec  and  v„  = 2650  m/sec. 

P2  ri2 


The  Kolymskaya  GES.  The  base  is  composed  of  quartz  plagioclase  granites 
(J-^  - Cr^ ).  The  average  thickness  of  the  belt  of  weathering  is  5 meters.  The 
unconsolidated  deposits  consisted  of  sand,  gravel  and  pebbles,  and  also  sandy 
loams  and  loams  (K  = 0 - 30  m).  The  thickness  of  the  permafrozen  rocks  was 

from  several  tens  of  meters  to  200-300  meters;  t = -6  to  -8°C.  The  mean 
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Table  1 Rayleigh  wave  velocities  in  the  region  of  the  Vilyuyskaya  GES,  the 
Kolymskaya  GES,  the  Deputatskiy  Mining  and  Concentration  Combine 
and  the  basin  of  the  river  Tatta 
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stones 
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values  are  v = 5^00  m/sec  and  vR  =3.15  km/sec.  The  section  of  the  housing 
P2  2 

settlement  is  composed  of  metamorphized  shales  (T  - J ) with  a northwestern 
course.  Their  belt  of  weathering  varies  within  considerable  limits  (from  5 
to  30  meters).  The  unconsolidated  deposits  (H  = 5-20  meters)  consist  of 
sand-gravel-pebble,  more  rarely  loams  and  sandy  loams.  The  seasonally  thawed 
layer  (in  the  working  period  up  to  2-3  meters)  is  mainly  loams  and  sandy 
loams.  The  basic  temperature  range  is  t^  = -1  to  -3°C.  The  mean  values  of 
the  elastic  wave  velocities  are  Vp^  = 5450  m/sec  and  vp^  = 2900  m/sec  along 
the  schistosity  and  v = 4600  m/sec  and  vR  = 2700  m/sec  across  the  schis- 
tosity.  P2  *2 

The  Deputatskiy  Mining  and  Concentration  Combine.  The  seismic  prospecting 
was  done  in  the  valley  of  the  Irgichan  River  (the  Indigirka  basin).  The 
bedrocks  are  Upper  Jurassic  metamorphized  sandstones  with  layers  of  clayey 
shales  and  siltstones.  The  belt  of  weathering  in  the  flood-plain  of  the 
Irgichan  River  has  not  been  established,  but  on  the  riverside  slopes  it  does 
not  exceed  2-3  meters.  The  unconsolidated  deposits  consist  of  sand-granite- 
pebble  deposits  and  also  of  sandy  loams  and  loams.  The  thickness  of  the 
permafrozen  rocks  in  the  region  of  the  work  is  200-300  meters  and  t = -6  to 
-8°C.  The  mean  values  are  v_  = 1+950  m/sec  and  v_.  = 2650  m/sec. 

p2  r2 

The  basin  of  the  river  Tatta  (the  Lena-Amga  interfluve).  The  bedrocks  are 
Middle  Jurassic  sandstones.  The  unconsolidated  deposits  of  the  section  of 
investigations  (H  = 40-45  meters)  are  composed  mainly  of  loams.  The  thick- 
ness of  the  permafrozen  rocks  exceeds  200  meters;  tp  = -4  to  -6°C.  The  bed- 
rocks are  characterized  by  a relatively  low  v^  of  *3800  m/sec,  which  can  be 
explained  by  the  higher  content  of  clay  particles  and  unfrozen  water  in 
the  sandstones.  The  values  of  vR  have  not  been  determined,  but  if  one  starts 
from  the  possible  values  of ^ 2 of  frozen  icy  sandstones,  for  such  a value 

of  v the  value  of  vR  x 2000^m/ sec . 

P2  2 
Discussion  of  the  Experiments 

1.  Comparison  of  the  above-presented  experimental  data  with  the  correspond- 
ing theoretical  dispersion  curve  and  the  dispersion  curve  for  the  data  of 
ultrasonic  modeling  shows  (Figure  8b)  that  the  theoretical  and  experimental 
curves  differ  at  H/\  = 0.1  - 0.8.  In  that  case  the  data  of  the  field  seismic 
survey  on  the  Rayleigh  wave  dispersion  (Table  l)  agree  well  with  the  results 
P*  of  three-dimensional  ultrasonic  modeling.  The  experiment  gives  a smoother 

drop  of  vR  from  the  value  v^  to  vR  with  increase  of  the  thickness  H of  the 
layer  on  the  half- space  2 1 than  follows  from  the  theory.  At  H/,\  > 

7 0.9  the  theoretical  and  experimental  curves  coincide.  The  greatest  scat- 
tering is  given  by  points  relating  to  the  base  of  the  Vilyuyskaya  GES,  and 
that  is  explained,  in  our  view,  by  the  relatively  small  length  of  the  survey 
lines  on  which  the  mean  values  of  v were  determined.  In  spite  of  the  fact 
that  bedrocks  of  the  above-considerea  regions  have  a different  petrographic 
composition  and  structure  and  are  characterized  by  different  (but  relatively 
* high)  values  of  v and  vR  (the  mean  values  are  v_  = 4. 6-5. 6 km/sec  and 

*2  2 2 
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Vp  = 2.65-3-2  km/sec)  the  dependence  v^/v^  for  all  those  rocks  is  a cocmcn 
2 * 2 
one,  which  makes  it  possible  to  recommend  its  use  in  the  practice  of  engin- 
eering seismic  prospecting  in  new  region.  These  are  above  all  regions  of  the 
propagation  of  igneous  and  metamorphous  rocks,  and  also  sedimentary  rocks, 
v and  v , in  which  they  fall  in  the  above  indicated  range . For  the  basin 
P2  R2 

of  the  river  Tatta,  where  the  sandstones  were  characterized  by  low  values  of 
v and  Vn  , vp/vn  = 0.75  at  H/\  = 1.25  agrees  with  the  theoretical  curve 
P2  2 U ri2 

for  the  corresponding  values  va  --  1.5  km/sec  and  v = 2.0  km/sec.  There- 

R1  p2 


fore  to  characterize  the  dispersion  of  v in  regions  of  the  propagation  in 
relation  to  weak  rocks  (sandstones  with  clayey  cement,  argillites,  etc) 
having  relatively  low  values  of  v and  Vp  in  a natural  bedding  (for  example, 
sandstones  of  the  Tatta  River  2 basin)  it  is  advisable  to  use 

the  calculating  dispersion  curve  obtained  for  specific  values  of  vQ  and  Vp  . 

It  should  also  be  noted  that  the  case  of  low-velocity  propagation  *2  ~2 

is  rare  for  regions  in  which  frost  is  widespread  and  occurs  when  in  the  above- 
indicated  rocks  water  and  air,  and  not  ice,  predominate  as  the  filler  of  cracks 
and  pores. 


2.  In  the  ultrasonic  modeling  an  investigation  was  made  of  the  influence  of 
the  loading  of  the  half-space  on  the  propagation  of  surface  waves  along  the 
boundary  of  the  layer  — the  half-space  (Figure  1 b,d,h,i).  The  dispersion 
of  the  surface  wave  velocity  was  studied  (Figure  8c).  In  the  presence  of  a 
layer  of  water  it  is  manifested,  starting  with  H/A  > 0.7,  and  Vp  is  reduced 
by  approximately  10%.  This  must  be  borne  in  mind  in  studying  Kthe  elastic 
properties  of  water- saturated  samples  of  rocks.  In  the  measurement  of  Vp 
in  them  it  is  necessary  that  the  water  layer  on  it  be  thin  (H  0.7 A), 

otherwise  we  obtain  reduced  values  of  v . And  to  exclude  in  addition  an 
intensive  wave  --  noise  on  the  water,  the  water  layer  above  the  sample  should 
be  reduced  to  a minimum  (less  than  1 mm). 


When  the  half-space  is  loaded  with  a layer  of  water-saturated  sand  and  ice 
a stronger  reduction  of  Vp  is  revealed  than  in  the  case  of  a layer  of  water 
(Figure  8c).  3y  comparing  the  obtained  dispersion  curves  with  one  another 
it  can  be  concluded  that  the  surface  wave  velocity  along  the  layer-half-spece 
contact  differs  all  the  more  from  Vp  (in  the  half-space)  the  greater  the 
shear  modulus  of  the  layer  (for  2 water  G =0  kg(force)/cm2,  for  water- 
saturated  sand  800  kg(force)/cm2  and  for  ice  ^ 29}0C0  kg(force )/cm2. 

3.  Comparison  of  the  seismograms  obtained  under  permafrost  conditions  [1,2, 
etc]  with  the  data  of  three-dimensional  ultrasonic  modeling  convinces  one 
that  both  in  absolute  values  of  the  longitudinal  and  Rayleigh  wave  velocities, 
distinctive  features  of  the  travel  time  curves  of  forward  and  refracted  waves, 
and  in  the  correlation  of  the  amplitudes  of  P and  R waves  and  in  the  correla- 
tion of  their  frequencies  a complex  picture  is  observed.  Complete  identity 
is  observed  in  the  correlation  of  the  kinematic  characteristics . With  respect 
to  dynamic  characteristics  of  the  wave  picture  it  should  be  noted  that  the 


14 


oscillograms  obtained  for  the  case  of  "ice  on  a half-space"  are  more  similar 
to  seismograms  for  the  case  of  "frozen  unconsolidated  deposits  on  a frozen 
half-space"  than  oscillograms  of  ultrasonic  modeling,  when  the  layer  was  of 
frozen  sand  or  frozen  loam.  For  the  latter  a strong  damping  of  longitudinal 
and  surface  waves  with  distance  was  observed,  whereas  recordings  on  an  ice 
layer  were  more  stable.  A reason  for  anomalous  lamping  in  sand  and  in  loam 
is  heterogeneities  of  the  structure  of  those  layers,  which  at  ultrasonic  fre- 
quencies can  be  commensurable  with  A.  The  more  homogeneous  structure  of  ice 
corresponds  better  in  that  respect  to  the  case  of  field  seismic  prospecting 
investigations  (the  size  of  the  heterogeneities  in  unconsolidated  deposits 
is  much  smaller  than  A).  The  data  of  modeling  and  seismic  prospecting  testify 
to  the  advisability  of  using  for  determination  of  S , G and /i  under  perma- 
frost conditions  Rayleigh  waves  characterized  by  ve^y  g£eat  intensity  on  the 
recordings.  High-frequency  recordings  of  R waves  can  be  used  to  calculate 
the  elastic  constants  of  frozen  Quaternary  deposits.  And  to  determine  the 
elastic  properties  of  rocks  hidden  under  alluvia  the  low-frequency  oscilla- 
tions should  be  registered  if  possible,  as  in  that  case  the  registered  value 
of  v^  differs  little  from  the  v^  of  the  rocks. 

4.  The  wave  picture  in  modeling  for  talik  conditions  was  compared  with  the 
corresponding  description  of  the  principal  types  of  waves  registered  during 
engineering  seismic  prospecting  [1,2].  Complete  similarity  was  established 

' with  respect  to  the  absolute  values  of  the  velocity,  characteristics  of  the 

time  travel  curves  and  correlations  of  the  amplitudes  and  frequencies  of  dif- 
ferent waves.  Under  conditions  of  thawed  rocks  the  main  information  about 
the  velocities  of  transverse  waves  in  rocks  hidden  under  a layer  of  alluvia 
or  water  is  borne  by  intensive  exchange  refracted  waves,  and  at  low  values 
of  H/A,  also  by  Rayleigh  waves. 

5.  In  three-dimensional  ultrasonic  modeling  in  the  presence  of  either  a 

solid  or  a liquid  layer  on  the  half-space,  in  subsequent  arrivals  a wave  with 
a boundary  velocity  v (the  velocity  in  the  half-space)  and  with  an  ampli- 
tude 2-k  times  as  ^2  large  as  that  of  the  P-j^Pi  was  registered.  Ac- 
cording to  the  conducted  investigation  it  is  a multiple  reflected- 

refracted  wave  which  has  been  reflected  from  the  boundary  of  the  layer  — 
air,  and  not  as  is  indicated-  hypothetically  in  [1]. 

6.  The  work  done  makes  it  possible  to  determine  more  precisely  a character- 
istic important  for  engineering  seismic  prospecting  — the  "depth  of  trapping" 

i_*.  . P of  elastic  oscillations  [8]  in  relation  to  transverse  and  surface  waves  in 

the  layer  on  the  half-space.  Actually,  the  value  of  H/A  on  the  dispersion 
curve,  starting  with  a certain  value  of  v /vR  , does  not  depend  on  H/A  , and 
that  means  that,  starting  with  that  valueK  2 of  H,  an  elastic  wave  does  not 
capture  the  half-space.  In  a recent  work  [8]  the  depth  of  capture  H„  is  es- 
timated by  A.  I.  Savich  at  0.25\  irrespective  of  the  type  of  waves  and  the 
elastic  characteristics  of  the  studied  media.  Analysis  of  the  theoretical 
and  experimental  dispersion  curves  obtained  by  us  makes  it  possible  to  con- 
clude regarding  H_  for  R and  S waves  in  a layer  on  a half-space: 
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l)  if  v_  (in  the  half-soace)  is  larger  than  vc  (in  the  layer)  by  an  crier 

2 , 1 
of  magnitude  or  more,  then  H,  ~ (0.2  - 0.2) X , which  coincides  -with  the  above- 

indicated  value  0.25a;  2)  if  v /v  = 5 - 10,  then  H-,  ~ (0.3  - 0.5)^  ; 3) 

"2  S1  J 

if  v /v  = 3 - 5,  then  H,  m (0.5  - 0.8)  A ; *+)  if  v /v„  = 1.1  - 2,  then 

S2  31  i S2  1 
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A. 


7.  During  ore  and  petroleum  seismic  prospecting  in  a permafrost  region  sur- 
face waves  usually  are  regarded  as  noises  and  are  not  used  in  the  interpre- 
tation. It  also  is  known  that  the  thickness  of  unconsolidated  deposits  under 
permafrost  conditions  is  determined  with  great  errors  in  many  cases.  Then, 
in  essence,  the  question  arose  of  determining  the  thickness  of  the  frozen 
layer  on  a frozen  half -space.  The  use  of  information  which  can  give  disper- 
sions of  the  Rayleigh  wave  velocity  to  solve  that  problem  ’will  make  it  pos- 
sible to  construct  seismogeological  profiles  more  substantiated^.  For  that 
purpose  we  will  use  a dispersion  curve  (Figure  8b)  characteristic  of  most 
regions  of  permafrost.  If  we  know  the  mean  value  of  v^  for  the  rocks  of  the 
given  region,  on  the  basis  of  values  of  v^  and  A-  1 2 determined  during 
seismic  prospecting  and  the  calculated  value  of  v /v  > we  find  the  value 
of  H/A  and  later  of  H.  The  value  of  H is  the  mean  1 R2  for  the  interval 
of  the'  seismic  prospecting  profile  on  which  v^  was  determined.  The  length 
of  such  an  interval  must  be  at  least  .100  meters.  In  order,  in  determining  H, 
to  deal  with  the  region  of  the  dispersion  curve  with  H/A  = 0-1,  Rayleigh 
waves  with  a wave  length  Ar,^  H should  .be  registered,  where  H_„„  is  the 

• ft  nicLX  -HcLX 

maximally  possible  thickness  of  the  frozen  unconsolidated  deposits  on  the 
section  of  the  work. 


Principal  Conclusions  and  Recommendations 

1.  In  determining  the  elastic  properties  (E  , G and  u ) of  unconsolidated 
and  solid  rocks  in  permafrost  regions  by  seiroiic^prospecting  it  is  advisable 
to  use,  besides  longitudinal,  Rayleigh  waves  (on  frost)  and  exchange  waves  (in 
taliks).  In  the  presence  of  a layer  of  frozen  unconsolidated  deposits  on  bed- 
rocks it  is  necessary  to  take  into  consideration  the  dispersion  v as  a func- 
tion of  the  layer  thickness  and  the  length  of  the  used  waves.  Understatement 
of  the  dispersion  leads  to  greatly  reduced  values  of  E^  and  G^  and  overstated 
values  of  u of  the  rocks. 

CL 

2.  To  determine  the  elastic  constants  of  the  massive  rock  covered  by  frozen 
alluvia  one  should,  if  possible  change  to  the  registration  of  low-frequency 
Rayleigh  waves,  since  in  that  case  v^  differs  little  from  v^  (in  the  rocks). 

2 

3.  In  the  investigation  of  various  questions  about  the  application  of  seismic 
prospecting  in  a permafrost  region  the  ultrasonic  modeling  of  wave  fields  on 
three-dimensional  models  of  natural  materials  (granite,  ice,  sand,  loam,  etc) 
is  promising.  Good  agreement  has  been  obtained  of  kinematic  and  dynamic 
features  of  the  wave  picture  according  to  the  data  of  modeling  and  field 
seismic  prospecting  work. 
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k.  The  Rayleigh  wave  in  a layer  on  a half-space,  the  dispersion  v„  of  which 
is  used  to  calculate  the  elastic  constants  of  rocks,  under  frost  conditions 
has  an  amplitude  5-20  times  as  great  as  that  of  a longitudinal  wave  only  at 
a layer  thickness  H <■  (0 . 3 - 0.5)A. 

When  the  layer  thickness  is  greater  A~  is  commensurable  with  and  the  phases 
of  a Rayleigh  wave  with  a large  amplitude  (10-20  times  greater* than  A ) bears 
information  about  the  Rayleigh  wave  velocity  in  the  layer.  In  the  presence 
of  a layer  with  H ^ A/3  close  to  the  radiator  only  a Rayleigh  wave  is  regis- 
tered in  the  layer,  and  a dispersion  of  the  phase  velocity  is  observed  at  a 
certain  distance  from  the  point  where  the  oscillations  are  excited  (for  con- 
ditions of  engineering  seismic  prospecting,  of  the  order  of  200  meters  or 
more ) . 


5.  In  the  presence  of  a thawed  layer  of  small  thickness  on  frost  an  intensive 
wave,  registered  in  "subsequent  arrivals,"  is  an  interference  (R  + 

wave.  In  the  presence  of  a layer  of  water  on  a frozen  medium  or  thawed 
naif-space  information  about  v„  is  born  by  intensive  waves  of  the  type  P^S^P-,  • 
Under  the  conditions  of  "a  thawed  layer  on  a thawed  half-space"  infor- 
mation about  vg  of  the  half-space  is  borne  by  P-^Sg?^  and  P^S 2S ^ waves  [7,3]. 

6.  Dispersion  curves  of  the  Rayleigh  wave  velocities,  calculated  with  for- 
mulas of  elasticity  theory  for  the  case  of  "a  layer  on  a half-space,"  give 

a more  rapid  drop  of  v^  from  v^  to  vR  at  H/A  = 0.1  - 0.8  than  follows  from 

the  data  of  field  investigations  in  the  rock  mass  and  ultrasonic  modeling. 

At  H/A  > 0.9  the  calculating  and  experimental  curves  coincide.  The  data  of 
seismic  prospecting  and  three-dimensional  ultrasonic  modeling  coincide. 


7.  The  dispersion  curve  vR/v^  (H/A)  obtained  from  modeling  and  confirmed  by 
seismic  prospecting  data  K2  is  recommended  for  determination  of  Vp  (in 
rocks)  in  most  permafrost  regions  for  which  the  mean  values  are  *2 

v > U.6  km/sec  and  v y 2.65  km/sec.  In  the  case  of  a low-velocity  half- 


space a calculating  dispersion  cur/e  for  specific  mean  values  of  v^  and  v^ 
of  the  given  region  should  be  used.  1 2 


8.  In  permafrost  regions  the  transition  from  a section  of  frozen  rocks  to  a 
talik  is  accompanied  by  a sharp  reduction  of  v . The  drop  of  vR  is  stronger 
the  greater  the  thickness  of  the  covering  bed  SH. 


9.  An  approximate  estimate  of  the  value  of  E of  permafrost  icy  rocks  can  be 
made  on  the  basis  of  measurements  of  Vp^  in  ^two  ways: 

a)  having  taken  mean  values  of  u^  equal  to:  0.3  (igneous  and  metamorphous 

rocks).  0.20  (quartz  sandstones),  0.26  (arkose  sandstones)  and  0.35  (lime- 
stones); 

b)  on  the  basis  of  charts  of  the  connection  of  E and  v of  rocks  of  a given 
region  (the  values  of  vp  and  E necessary  for  construction  of  the  diagram  are 
determined  during  seismoacoustxc  work  in  adits  and  wells  and  on  outcrops  where 
the  registration  of  vg  or  v^  presents  no  difficulties) . 
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10.  During  measurement  of  v_^  in  water- saturated  samples  of  rocks  it  is  ne- 
cessary that  the  water  layer*lon  them  be  minimal  (better  up  to  1 mm)  for  reduc- 
tion of  the  -dispersion  v^  and  that  recordings  not  be  complicated  by  an  inten- 
sive wave  in  water. 

11.  It  is  advisable  to  use  the  dispersion  of  Rayleigh  waves  in  a permafrost 
region  also  in  the  construction  of  seismogeological  profiles. 
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